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ABSTRACT 


The objective of this work has been to enhance the predictive capability of widely used CFD codes 
through the use of solution adaptive gridding. Most problems of engineering interest involve multi- 
block grids and widely disparate length scales. Hence, it is desirable that the adaptive grid feature 
detection algorithm be developed to recognize flow structures of different type as well as differing 
intensity, and adequately address scaling and normalization across blocks. In order to study the ac- 
curacy and efficiency improvements due to the grid adaptation, it is necessary to quantify grid size 
and distribution requirements as well as computational times of non-adapted solutions. Flowfields 
about launch vehicles of practical interest often involve supersonic freestream conditions at angle 
of attack exhibiting large scale separated vortical flow, vortex-vortex and vortex-surface interac- 
tions, separated shear layers and multiple shocks of different intensity. In this work a weight func- 
tion and an associated mesh redistribution procedure is presented which detects and resolves these 
features without user intervention. Particular emphasis has been placed upon accurate resolution 
of expansion regions and boundary layers. 

Flow past a wedge at Mach = 2.0 is used to illustrate the enhanced detection capabilities of this newly 
developed weight function. Figure 1 presents weight functions evaluated using the previous proce- 
dure, lower half plane, as well as the current procedure, upper half plane. 
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Figure 1. Comparison of Weight Functions. 


ADAPTED GRID WlTH/WlTHOUT (BOOLEAN) 
(M.50.5) 



Figure 2. Comparison of Adapted Grids. 


It can be observed that both weight functions clearly detected the primary shock. It can also be seen 
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that the expansion fan, boundary layer, and the reflected shocks are much more clearly represented 
in the current weight function. Adapted grids using both weight function formulations are presented 
in Fig. 2. The high gradient regions of the expansion region are only reflected in the adapted grid 
using the new weight function. The reflected shock is also much sharper. Figure 3 compares the 
solution obtained using the current adaption procedure with that obtained using the original grid. 
The enhanced resolution is clearly evident. 


ADAPTED/ORIGINAL SOLUTION (BOOLEAN) 



IJO> MJO UW 

Figure 3. Comparison of Solutions Using Adapted Grid. 


Supersonic flow at Mach= 1 .45 and 1 4 degree angle of attack has been simulated around a tangent- 
ogive cylinder. The grid and associated flow solution constructed after two adaption cycles using 
hybrid differencing of the grid equations and the current weight functions is presented is presented 
in Figure 4. 



Figure 4. Adapted grid after two cycles. Figure 5. Adapted grid after two cycles. 


Figure 5 presents the grid constructed using the previous weight function and the same flow condi- 
tions and number of adaptation cycles. Figures 5 and 6 present streamwise cuts of the two grids 
shown in Figs 4 and 5 at X/D = 5.5 and 7.5 respectively 
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Figure 6. X/D = 5.5 


Figure 7. X/D = 8.5. 


Figure 8 present the flow solution obtained using the NPARC [NASA 1993] flow solver, the KE 
turbulence model option and two adaptation cycles. Figure 9 presents the associated weight func- 
tion. 



Figure 8 Normalized Stagnation Pressure. 



Figure 9. Weight Function. 


Examples will presented to demonstrate the capability for solution-adaptive regridding of multi- 
block launch vehicle simulations. 
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COMPLEX GEOMETRY / COMPLEX PHYSICS 
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4. Robust redistribution procedure for use with weight function. 


GOVERNING EQUATIONS FOR GRID MOVEMENT 
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EVALUATION OF FORCING FUNCTIONS 
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P k =W§ k /W, k= 1,2,3 
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• Minimum variation of coefficients. 



EVALUATION OF WEIGHT FUNCTIONS 
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maximum in order to maintain the relative strength. 


OVERALL SOLUTION PROCEDURE 
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SOLUTION OF GRID EQUATIONS 
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Adapted Grid Central Di 



1052 



H 

£ 

W 


o 


H 

Z 


00 

PQ 

Pi . 

§ § 
^ 0-> 
P o 

G cd 


op 5 

^ys 

PQ Vh 

__, C3 

cd £ 
Ch *g 

o -P 
T* c3 

/5 ^ 

PP o 

g a 
S o 

rP cd 


S ca "t2 *-> 

S C/3 * ^ f-H 

O <D C3 7? 

W) es >=t <D 

■S 5 ^ 

o £ o cL 

• • Z 2 


3 « 

rS o 

C ,3 

+-> ^1— I 

00 *-h 

« ° 

o c 

43 .2 

a s> 

M £2 

>. T3 
b <D *J 

V»NJ • 1-H * ^ 

<4—1 O 
22 -»-h ’23 

g o a, 

3 <D * 

O P, W 

PP ^ • 


1054 


• Local iteration tor desired distribution. 

• Can be used to keep sharp comers, and to transfer information 
between blocks. 
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7 KE Model 




Case 1 NP ARC Grid 7 a2 KE Model 
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Improved convergence behaviour. 

More closely resembles experimental data. 



ONGOING WORK 
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y parabolized gnd equa- 



ADAPTED DISTRIBUTION MESH 
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Scramjet Inlet NPARC Solution (lower), Weight Func (upper) 
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